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Abstract 

We examine again the problem of the damping rate of a moving heavy fermion in a 
hot plasma within the resummed perturbative theory of Pisarski and Braaten. The ansatz 
for its evaluation which relates it to the imaginary part of the fermion propagator pole in 
the framework of a self-consistent approach is critically analyzed. As already pointed out 
by various authors, the only way to define the rate is through additional implementation 
of magnetic screening. We show in detail how the ansatz works in this case and where 
we disagree with other authors. We conclude that the self-consistent approach is not 
satisfactory. 
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1. Introduction 



Hot gauge theories are the framework to study perturbatively plasmas of weakly in- 
teracting particles. Of special interest are the properties of quasiparticle modes. The 
knowledge of the fermion self-energy, first computed by Klimov and Weldon [1] , has yielded 
the dispersion relations of thermal fermionic excitations at leading order in the coupling 
constant g (in hot QED/QCD). They are obtained by equating the real part of the inverse 
effective propagator to zero. 

To go further in the quasiparticle interpretation, a question has arisen : what is the 
damping rate of this quasiparticle ? In other words, by which coefficient T the large time 
e~ rt behaviour of the corresponding response function is governed ? When the effective 
fermion propagator has a pole, an hypothesis always assumed in previous works [2-7], the 
width T of the quasiparticle is given by the opposite of the pole imaginary part. 

Initially, the damping rate of a gluon at rest raised a problem : one-loop calculations 
led to gauge-dependent results. The solution of this "plasmon puzzle" required to take 
into account all diagrams contributing to leading order in g [3, 8], thus giving a finite and 
gauge-independent result [8, 9]. 

Similarly, a finite damping rate V was found for a massive fermion at rest [2]. In 
that case, only the longitudinal component of the gauge boson propagator contributes to 
r and the long range interactions are screened by the Debye length in the framework of 
the resummed theory [10]. As soon as the fermion is moving, the transverse component 
also contributes. Its behaviour in the space-like region leads to a logarithmic infrared 
divergence [2]. The resummation program is apparently not sufficient to screen the long 
range magnetic interactions. 

Many authors [3-5, 7] have nevertheless tried to calculate V with the hypothesis that 
the fermion propagator has a pole, i.e. to find the pole imaginary part. Their attempts to 
screen the infrared divergence are based on a self-consistent approach, originally proposed 
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by Lebedev and Smilga [3], but results are often different or even contradictory. Our 
purpose here is to study this situation. In agreement with previous authors [4], we find 
that pursuing literally the search of a pole of the propagator (i.e. naively using a "complex 
energy-shell" condition to define the rate) in the framework of the self-consistent approach 
leads to a failure : the logarithmic singularity is not screened. On the other hand, going on 
with this approach in the QCD case, in presence of magnetic screening, leads to an answer 
for the rate only in the case of a fast moving fermion ; the corresponding equation for a non 
relativistic heavy fermion has no solution. The result obtained for the relativistic fermion 
differs from the usual intuitive "narrow width" expression, 0(g 2 T£n(l/g)) pointing out at 
possible difficulties connected with the simple minded analytic continuation performed to 
define the propagator in the complex energy plane. In this respect, let us stress that our 
results disagree with those of ref. [7] . We think that the source of the disagreement lies in 
the way analytic continuation is performed to evaluate the damping rate. 

After relating the damping rate with the presumed pole of the fermion propagator 
(section 2), we calculate the fermion self-energy in the imaginary time formalism and 
perform the analytic continuation to complex external energy (section 3). In section 4, 
first discussing the solution obtained without magnetic screening, we then analyze the QCD 
case in presence of a magnetic mass both for a nonrelativistic and a fast heavy fermion. 
Section 5 is the conclusion. 

2. Definition of the damping rate 

We assume the existence of a pole in the complete fermion propagator. The damping 
rate is then given by the opposite of the imaginary part of this pole. Working in the 
imaginary time formalism, the complete fermion propagator is : 

S(P4,P) = . T - „ 1 „ - ( ^ , (1) 
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where T, E is the fermion self-energy in euclidean space*. The usual analytic continuation 
for the retarded propagator is ip^ — > po + is with p real. In order to locate the position 
of the complex pole in the lower half energy plane, we tentatively assume that we may 
continue the retarded propagator to complex values of po as ip4 — > po- The pole is the 
value of po for which : 

det{po7o-Pl-M-Z{p ,p)} = , (2) 
where H(po,p) = — Ee(— zp o ,£0- In general, po is complex : 

Po = E-iT , (3) 

with T the fermion damping rate. 

By rotational invariance in the rest frame of the plasma, E takes the form : 

E = a7o + bp 7 + c , (4) 

which from eq. (2) leads to : 

po = a± v / (l + b) 2 p 2 + (M + c) 2 . (5) 

In eq. (5), a, b and c are functions of po so that eq. (5) is a transcendental equation, 
yielding the pole po = E — iT as an implicit function of \p\ . 

For a heavy fermion (M >> T), a « y/p 2 + M 2 and Re p = E ~ a/p 2 + M 2 . We 
use the perturbative approximations 6 << 1, c << M, and neglect quadratic terms in b 
and c. Then we show that : 

* In the imaginary time formalism, the euclidean Dirac algebra {7 M , 7^} = 26^ is used. 
We note = p^lo + P 7, where P4 is the discrete euclidean coordinate p^ = (2n + 1)ttT. 
The continuation for 7-matrics is achieved by — > —7 (70 unchanged) to recover the 
algebra {7^,7^} = 2g' iV . 
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-Im po ~ --^-Im Tr(J) J r M)T,(p ,p) 
AE 



(6) 

P o=E-iT 



where = E-f — p 7. Notice that T is determined by a functional equation since E is 
evaluated at po = E — iY. In principle, going off the real energy axis requires the knowledge 
of the analytic structure of E, namely the existence and positions of cuts if any, as well as 
the different sheets, especially the one on which the pole is to be searched. 

Since this analytic structure is largely unknown, it may be dangerous to explore the 
complex plane naively without any guidance, as will be discussed in sect. 3. 

3. Calculation of Im E for a complex external energy 

The leading contribution to the hard fermion self-energy in the imaginary time for- 
malism reads : 



Z(p) = -g 2 C F [ ^\ T Mtfhu'^iq)} (7) 

J "f* (27r) q4 ^nT 



Cp is the QCD Casimir which should be replaced by 1 in the case of QED. Sq(p') is the 
bare internal fermion propagator (p' = p — q) : 



S oiP') = tL 4^ with £' 2 = ? 2 + M 2 (8) 



*A^ u (q) is the resummed gauge boson propagator including hard thermal loops [11]. In 
the Coulomb gauge, it consists of two terms : 

• the longitudinal part *A 00 = *A e (q). 

• the transverse part *A lj = Q lj *A t (q), where Q lj = 5 lj —q l q j , with % = j^. 
*Ag gives a finite contribution to Im E, whereas *A t is responsible for infrared problems 
[2] . Therefore we shall focus on the contribution which comes from the transverse part : 
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Z t = -g 2 C F [ ^Tj2*MQ±,Q)WSo( P 'h j Q ij } ■ (9) 

J soft {^r x 1 J 

In order to sum over q^, we use the spectral density representation for the transverse 
propagator : 

rP r°° 

*A t (q 4 ,q) = dre^ T du p t (oj, q) (1 + n(u)) e~ u T , (10) 

JO J-oo 

where n denotes the Bose-Einstein statistical factor and the transverse spectral density 
Pt(co, q) given in [11] is proportional to the imaginary part of *A t ((74, q) after the continu- 
ation i(/ 4 — > u + is : 

p t (w,q) = -Im*A t (q 4 -> -i(u) + ie),q) . (11) 

7T 

To calculate Im E for a complex external energy we may use a mixed representation for 
the bare fermion propagator : 

S (p) = dr e ipiT [n(-E) e- ET P+ + h(E) e Er P_\ (12) 
Jo 

h refers to the Fermi-Dirac statistical factor and P + (P_) project on positive (negative) 
energy states : 

P ± = ^[E 7o ±(tP7 + M)] . (13) 

Instead of using eq. (12), Pisarski [7] proposes an ansatz for the fermion propagator in the 
spirit of eq. (10) : the fermion spectral density is taken as a Breit-Wigner distribution. 
This however does not allow exploring the pole position in the lower half complex plane. 
The analytic continuation which is used in [7] in order to go from the evaluation of Im £ 
on the real axis to its value on the complex pole seems quite questionable to us as we shall 
argue in section 4. 
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We insert eqs. (10) and (12) in eq. (9) and find after the proper analytic continuation 
(keeping only the soft space-like region \q \ < q << T) : 

Im E(p ,p) = 2g C F T / -—3 / — Pt ( qo ,q)I m — 2 — (14) 

J soft l /7 U J-q Q0 Po - & 

which is the result found in [4] with the real time formalism. We have used \qo\ << T << 
E' and n{qo) ~ ^. Moreover, as po is hard, \po\ ~ E, we replace p' = p — qo by E in the 
numerator of eq. (14). We also write : 

1 = lfJ_ + J_LlJ- 
p' 2 -E' 2 2p' \p' -E ' p' +E'J 2E p' -E> 

as E and E 1 are hard energies. We get : 

Im S( P0 ,p1 = 9 2 CS [ ^ / +<? ^p t (q , q) (e 1o - {p'Wl) ~ m) Im - 7 -^— 

^ J soft {^) J-q Q0 V J PO - E 

(16) 

From eqs. (6) and (16) we obtain : 

r (P) = /i 2 u / d 9 ? / Pt(qo,q) dx(l + -x - x )Im — 17) 

where u = £ is the fermion velocity and x = cos 9 with 6> the angle between q and p. The 
cut-off q* determining the soft integration region is some arbitrary scale «T e.g. Ty/g. 

4. Calculation of the damping rate 

In order to compute T(p), the "narrow width" approximation was used as a first 
attempt ; as in the case of a fermion at rest, the external energy was taken to be real. But 
this led to an infrared logarithmic divergence [2] . In order to cure this pathology, Lebedev 
and Smilga [3] proposed a self-consistent approach according to which the non-zero value 
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(15) 



of r itself would screen this divergence. The practical implementations of this philosophy 
differ however in the literature. Here we first follow the procedure of Baier, Nakkagawa 
and Niegawa [4] which allows us to deal with real or complex external energies. We now 
restrict ourselves to p > T and neglect the % term in eq. (17). 

The self-consistent approach amounts to introduce the rate V as a dissipative coeffi- 
cient in the bare fermion propagator [12], so that in eq. (8), E' is replaced by E' : 

E' = E'-iT{p') . (18) 

In other words, the only singularity of the effective dissipative propagator is the complex 
pole at E — iT. With these modifications, eq. (17) becomes : 

g 2 C F T 2 /■«* 2 H dq 
F(p) = . 2 v / dqq Pt{qo,q) J{qo,q) 

47r J0 J -a QO 



with 



J( c lOj ( l) = [ dx(l — x 2 )Im [ — 

j-i yp'o-E'^ 



(19) 



Approximating T(p') ~ T(p), one finds [4] 



g 2 C F T n* H dq 

r (p) = o 2 v dqq pt{qo,q) 

^ Jo J- q Qo 



fRep -E + vq-q \ fRep -E-vq-q \\ 

arctan —— — arctan —— > . (20 

V Imp + r(p) J \ lmp + T{p) J J V ' 

To apprehend the infrared behaviour of the r.h.s. of eq. (20), it is sufficient to put q = 
in the curly bracket, and to use the limiting form of Pt(qo, q) in the space- like region [11] 

= 1 1 I 4 m 9 qi 

Pt{qo,q) — ; -g , 21 

q «q«m g q 2 ^ 1+ h^ m2 qo\ 

where m g is the thermal mass of the gauge boson. We get : 
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= - 5^ v [ — — arctan [ ^rn—^- I • 
U ) 2tt 2 J q TV \4 q 2 J 

( ± Re po - E + v q ± Re p - E - vq\ 

. < arctan— — — arctan— — . > 

\ lmp +T(p) lmp Q + T(p) J 

a)Absence of magnetic screening (QED case) 

As noticed in [4], with po = E — iT, which is required in principle by eq. (6), the r.h.s. 
of eq. (22) presents a singularity. Thus the self-consistent approach alone does not achieve 
the initial aim which was to screen the logarithmic divergence. In fact the origin of this 
divergence is well-known : it is due to the factor in eq. (21) which reflects the absence 
of magnetic screening in the static limit at this order. It may seem surprising to try to 
cure a bosonic disease by acting on the fermionic sector, as in the self-consistent approach, 
and the success of such a treatment in the screening of the divergence would have been 
somewhat miraculous. 

Let us stress that the results found in [5] and [7] for the non-relativistic fermion by 
putting Im po = in the r.h.s. of eq. (22) do not in principle give the damping rate, 
since this procedure deviates from the prescription of eq. (6). Although a dissipative 
coefficient is introduced, the choice Im p° = is a kind of reminiscense of the narrow 
width approximation, the validity of which is precisely questioned by the divergence at the 
pole : thus the link with the damping rate is unclear and the fmiteness of the resulting 
expressions may be misleading. 

Strictly speaking, the r.h.s. of eqs. (20), (22) rely on the replacement T(p') — > T(p) 
in eq. (19). Since V is not independent of p, this approximation is correct as long as 
\T(p') —T(p)\ « |Im p°+T(p)\, but it fails at Im p° = — T(p), so that the conclusion of [4] 
is not rigorously stated. Since T(p') — T(p) ~ O(q) the crude replacement of Im p + T(p) 
by this 0(q) term in eq. (22) still does not cure the divergence at the pole : the conclusions 
of [4] are expected to remain unaltered. 
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(22) 



We will show however that the situation worsens : we obtain a functional equation 
for T(p) which has no solution even in presence of magnetic screening (at least in the non 
relativistic case) ! Using 



dT 

r(p') ~ r(p) - q — x 



(23) 



which holds* if 
is given by : 



q dT 

r d P 



« 1, and with EE' = E + |^ - vqx, the function J in Eq. (19) 



I (l- X 2 )(l m p 0+ T(p)-q§x) 

dx — 

i [Re po - q - E - |J + vqx} 2 + [Im p + T(p) - q^ x] 2 



(24) 



Putting p Q = E — iT, and after some simple manipulations we find (in the perturbative 
approximation ^ << v) that J may be safely approximated by : 



J 



Then eq. (19) gives : 



go g 
q 2E 



(25) 



go _g_ 
g 2E 

(26) 



* Eq. (23) assumes the differentiability of T(p) with respect to p. This assumption may 
seem quite restrictive and not founded. Let us notice however that : 

1 - it is natural in the framework of the pole ansatz which supposes that the pole is 

isolated from other singularities. 

2 - it is less stringent than assuming T(p') = T(p). 

3 - Finally it is more consistent than the hypothesis of constant behaviour, which even- 

tually leads [5, 7] to T(p) being a differentiable function of v(p) = p/E(p). 
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As far as the QED case is concerned (absence of magnetic mass) the behaviour of pt(qo, q) 
in the space-like static limit (cf. eq. (21)) still generates a logarithmic singularity. 



b)Presence of magnetic screening (QCD case) 

The unexpected novelty is that the occurrence of magnetic screening does not seem 
to solve the fermion damping puzzle, at least in the non-relativistic regime ! 

Unfortunately we do not know the spectral function pt(qo, q) in presence of a magnetic 
mass, except in the static limit. In order to explore the effect of magnetic screening, we 
propose a simple parametrization which interpolates smoothly between : 

• the static limit, where we expect the magnetic "mass" p to screen the infrared 
singularity : Re U t (q ,q) p 2 so that pt{qo,q) 90 -o ni } ni , instead of X as in 
eq.222222221). 

• the region away from the static limit, where the hard- loop approximation of pt(qo, <?) 
certainly holds. Hence we take : 

3 m 2go 

P?° del (qo,q) = —-2 • (27) 



This leads to 



(q 2 + V 2 ) 2 + {h m D 2 (f 



= —sgn I — — J - F — v [ dq 



dp J An 2 J q 2 + p 2 

( TP" \ f( v ~ m) \ f( v + w) 

2 arctan 2bj . + arctan -^ — *2J-v - arctan -^ — 2J ^- 

\q 2 + p 2 J \ q 2 + p 2 I \ q 2 + p 2 



(28) 



where v = f^m^. 



One can show that the curly bracket in eq. (28) is dominated by the first term yielding: 



dT\ g 2 CpT f°° r ( r 

r(p) = -sgn ( —j ~^~v ^ dr-^arctan [K-^ ) . (29) 
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where we have set r = ^ and K = §7T^. The g* dependence of eq. (28) is subleading in 
g and can be dropped. 
The result is : 



r(p) = -s#n — 



\dp J An 



-v sink' 



3tt 

16 E/x 



(30) 



Consequently T(p) should satisfy the condition 



dr 

dp 



P 



g 2 C F T d 

47T dp 1 ^/ p 2 + M 2 



sink 1 



3tt 



(31) 



When the expression between brackets is an increasing function of p, eq. (31) has no solu- 
tion. This happens when p < M. Thus there is no solution for a non-relativistic fermion: 
the corresponding propagator has no pole. In the case where v « 1, the qo integration 
range is restricted to the static limit region q << q where the parametrization (27) is 
certainly reliable. When p > M, eq. (31) has a solution. Namely, in the ultrarelativistic 
case we get : 



T(p) \ v=1 = 9 CpT sink 1 

47T 



3tt m z g 
16 



(32) 



Assuming fi ~ g 2 T, ~ |j << 1, we thus have 



r(p) 



v=l 



3g 2 C F T m 2 g 
64 E/x 



(33) 



Surprisingly, we do not find the logarithmic dependence in \i naively expected since \i 
screens the logarithmic divergence. We may first question the reliability of the parametriza- 
tion (27) used in the whole go integration range. However, due to the sign flip of J(qo, q) 
and the odd property of Pt(qo, q) under the change qo —qo, the essential part of the qo 
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integral in eq. (26) comes from the narrow window [— ^] where the ansatz (27) is 
expected to hold *. Then the nature of the /(/-dependence of T(p) (for p > M) is controlled 
by the ratio of the width of the window (-|) to the width of the spectral density p™ odel 

2, 2 

considered as a function of m which is ( 1 M 2 )• The magnitude of this ratio R given by : 

1 4 7rm g 

R = K^— , (34) 
r 2 + 1 

with r = ^, is controlled by the parameter K defined in eq. (29). 

m 2 

If p is arbitrarily small, namely p « -j=r « 0(g 2 T) the spectral density pt is 
sharply peaked. It is then legitimate to replace as done in previous works [4], q by in 
the integrand of eq. (26), except in Pt(QoiQ)- In this case, instead of (33), eq. (32) 
yields indeed the logarithmic dependence on p : 

But on the other hand, for the physically sensible value p ~ 0(g 2 T), K « 1 and the 
density p t has a broader support than the window — ^ < ^ < This forbids the 
appearance of the logarithmic dependence and leads to eq. (33). 

We point out that the sign flip in J, hence the narrowness of the window instead 
of 0(1), are a direct consequence of the self-consistent approach associated with the pole 
prescription of eq. (6), through eqs. (23)-(24). If one would instead use a nondissipative 
fermion and work with a real energy p° - as usual in the narrow width approximation -, 
J would be replaced by the integral over cos 9 of a 5-function, which would not produce 
such a sign flip. In such a case the corresponding ratio R would be controlled by m 2 ,/// 2 , 
much larger than 1. One would then get the standard £n(m g /p) dependence in T. This 



The second window let open by J and the odd behaviour of p t , which is v — ^ < 
V < v + 2fJ, gives only a negligible contribution with respect to the one we discuss. 
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remark leads to the conclusion that even in presence of magnetic screening and when the 
pole exists, the narrow width approximation and the self-consistent approach disagree. 

This unconventional conclusion, associated with the qualitative difference between the 
nonrelativistic and ultrarelativistic regimes without any apparent physical reason make us 
feel uneasy. Hence we wonder about the correctness of the self-consistent procedure, i.e. 
the pole ansatz of the dissipative propagator together with the exploration of the lower 
half complex energy plane, irrespectively of the true analytic structure of the complete 
fermion propagator. 

In ref. [7] , arguments are given to obtain informations about the singularity structure 
of the propagator and discuss the possibility of ignoring this structure, keeping only the 
pole. We do not, however, think that these arguments are well founded. In fact, the 
analytic continuation which is used in ref. [7], in order to explore the lower half-plane 
of the energy, is performed on the self-energy imaginary part evaluated for real external 
energies. But, clearly : 



Im S(real po — > complex p ) ^ Im (E(real po — > complex p )) 

and it is only the latter quantity which is of interest and which we actually study in the 
present work. 

In fact, instead of eq. (12), ref. [7] uses a spectral density representation for the 
internal fermion propagator, which is correct only for a real external energy p . To take 
into account the width of the fermion, a Breit-Wigner form for the fermion spectral density 
is introduced by hand. To make contact with this method, we write the imaginary part of 
eq. (15) with a dissipative energy E' as : 

x 1 r(p')+Im po 

m p' — E' (E-E>- q o r + (IV) + Im p f [M> 

We recover the same Breit-Wigner form as ref. [7] only for Im p = 0. The ad hoc analytic 
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continuation which is proposed in order to go to Im p = —V leads to results contradicted 
by our treatment, which enables us to go directly to the pole. In particular, this is flagrant 
for the case of a non-relativistic fermion for which ref. [7] states that no magnetic screening 
is needed. Also when going to the complex pole in eq. (36), one finds as a consequence that 
the final answer for T{p) cannot depend on itself as appears in [7] but only on T(p') — T(p) 
which is what we have dealt with. 

5. Conclusion 

By a careful exploration of the location of the pole of the fermion propagator, we have 
shown that the self-consistent approach using a dissipative fermion is of marginal relevance 
to solve the heavy moving fermion damping puzzle : 

• it does not provide any screening for the infrared divergence, which has to be cured 
by magnetic screening. 

• In presence of magnetic screening, the situation depends on the kinematical regime. 
In the non relativistic regime the functional equation obtained for V has no solution. 

The situation changes dramatically in the relativistic regime, but without any apparent 
physical reason : we do find a pole in this case, hence a value for T. However we got a 
rather unusual result, especially concerning the dependence of V on the magnetic mass //, 
as compared to the commonly expected In fx. The /(/-dependence which we get in V is 
logarithmic only if /i is very small with respect to the standard scale 0(g 2 T) and this is 
a direct and specific consequence of the use of a dissipative fermion. This disagreement 
with the result obtained in the narrow width approximation (which is naively expected to 
work since F ~ g 2 T « gT << M) suggests that the self-consistent approach is not a safe 
procedure and may lead to question the validity of the simple hypothesis of a leading pole 
singularity. The complete singularity of the propagator is unknown : there are no physical 
requirements to determine it and in particular no way to define a physical sheet as for the 
T = S'-matrix element. 
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Notice that the absence of a pole in the fermion propagator does not jeopardize, in 
principle, the quasi-particle interpretation of plasma excitations. Indeed, a quasi-particle 
shows up in the large time asymptotic behaviour of the mixed time/momentum represen- 
tation of the (retarded) fermion Green's function as : 

G(t,p) tS+oo A(t,p) exp[-iu (p)t - F(p)t] 

where u)o(p) gives the dispersion relation and T(p) is the damping rate. In the complex 
energy plane, po = LVo(p)—iT(p) is the nearest singularity of the propagator with respect to 
the real axis, as can be shown, e.g. by a saddle point approximation, and this singularity 
need not be a pole, but e.g. a branch point as claimed by Smilga [13] in the QED case. 
The prefactor A is non exponential in time and depends on the nature of this singularity 
(for example it is a constant for a pole). In practice however, locating any singularity 
off the real energy axis (a pole, a branch point or an essential singularity) involves the 
knowledge of the analytic structure of the finite temperature propagator which is out of 
reach at present. 
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